Interfacial Magnetoelectric Coupling in Tri-component Superlattices 
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Using first-principles density functional theory, we investigate the interfacial magnetoelectric cou- 
pling in a tri-component superlattice composed of a ferromagnetic metal (FM), ferroelectric (FE), 
and normal metal (NM). Using Fe/FE/Pt as a model system, we show that a net and cumulative 
interfacial magnetization is induced in the FM metal near the FM/FE interface. A carefully anal- 
ysis of the magnetic moments in Fe reveals that the interfacial magnetization is a consequence of a 
complex interplay of interfacial charge transfer, chemical bonding, and spin dependent electrostatic 
screening. The last effect is linear in the FE polarization, is switchable upon its reversal, and yields 
a substantial interfacial magnetoelectric coupling. 



Materials with coupled magnetic and electric degrees of 
freedom are classified as magnetoelectrics. A strong mag- 
netoelectric coupling enables control of magnetism by the 
electric field or vice versa, and hence has a strong appeal 
for emerging device applications ,iiSi2i4 A linear coupling 
between these two degrees of freedom involves breaking 
both space and time reversal symmetries and is therefore 
of fundamental interest. Multiferroic oxides in which si- 
multaneous ferroelectric and magnetic orders exist form 
a class of promising magnetoelectric materials. However, 
many single phase multifcrroics have rather limited ap- 
plication either because of their low Curie temperature 
or weak coupling. 

Heterointerfaces have proven to be ideal for control- 
ling and manipulating electrical charges and spins in 
solid state devices. Recently, efforts have focused on 
the so called "interfacial magnetoelectricity" in which 
a magnetoelectric coupling arises at a ferromagnetic 
metal/dielectric (or ferroelectric) int er f ace £iiLLS£ i i 1 1 1 1 1 1 2 
A ferroelectric or a dielectric (upon the application of an 
electric field) can induce free charges near the interface 
within a ferromagnetic metal. These screening charges 
are spin dependent in the ferromagnet, and thus in turn 
yield an additional magnetization that exists only within 
a nanometer near the interface. In contrast to the in- 
trinsic coupling in multiferroic oxides, this effect can be 
viewed as an extrinsic magnetoelectric effect and may of- 
fer an alternative to a single phase multiferroic in provid- 
ing the robust room temperature magnetoelectric (ME) 
effect. 

In this paper we explore "tri-component" FM/FE/NM 
superlattices consisting of alternating layers of a ferro- 
magnetic metal, ferroelectric, and normal metal as de- 
picted in Fig. [TJ Unlike the symmetric heterostruc- 
ture^'^ ' 10 i 11 ' 12 the broken inversion symmetry in an 
asymmetric superlattice permits a leading magnetoelec- 
tric coupling in the free energy expansion that is linear 
in the electric polarization of the FE, 8 i.e., oc PM 2 , and 
hence can be electrically controlled. Moreover, since the 
induced magnetization only occurs at the FM/FE inter- 
face, the overall induced magnetization does not cancel 
as in the symmetric structure^ and can cumulate in a su- 
perlattice to achieve a large macroscopic magnetization. 

Using density functional theory (DFT), we demon- 
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FIG. I: Schematics of the "tri-component" superlattices 
which consist of alternatively stacked ferromagnetic metal 
(Fe), ferroelectric, and normal metal (Pt). 



strate this effect using Fe/FE/Pt superlattices, where for 
FE we use BaTi0 3 (BTO) and PbTi0 3 (PTO) as two ex- 
amples. PTO and BTO are two prototypical FEs with 
robust room temperature polarization. Their significant 
difference in polarization allows us to investigate the po- 
larization dependence of the ME coupling effect. A rather 
complex picture emerges as a result of a detailed elec- 
tronic structure analysis. The interface magnetization is 
controlled by the interplay of quite different physical ef- 
fects including electrostatic screening, contact potential 
difference, formation of the Fe-0 chemical bonds, metal 
induced gap states, and subtle chemical differences re- 
sulting from the polarization switching. The electrostatic 
effect is approximately linear in the FE polarization and 
switchable upon the reversal of the polarization, in agree- 
ment with the analytical models The overall magneto- 
electric coupling we find in these systems is remarkably 
large compared to many other existing composite struc- 
tures that have been reported so far. 13 i 14 

We use DFT within the local density approximation 
(LDA) as implemented in the VASP code and projector 
augmented-wave pseudopotentialsJ^ We apply a plane 
wave cut off energy of 600 eV and the 8x8x8 and 8x8x2 
k-point meshes for the Brillouin zone integration for the 
bulk and superlattice, respectively. Calculations of the 
lattice parameter in the tetragonal phase yield 3.86A, 
c/a = 1.04 for PTO and 3.95 A, c/a = 1.01 for BTO. 
The lattice constants are about 1.5% less than the ex- 
perimental values of 3.904 and 3.994A^ as is typical for 
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FIG. 2: The planar (thin gray curve) and macroscopi- 
cally averaged potential (thick green curve) as a function 
of distance along the [001] direction for the fully relaxed 
(Fe)ii(PbTi0 3 ) 9 (Pt)8 superlattice. The straight line in PTO 
marks the residual field in the PTO unit cells, where an mag- 
nitude of ±0.06 eV/A is calculated. The arrows mark the 
direction of the FE polarization in PbTiOa. 
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FIG. 3: The rumpling parameter for each PbTi03 atomic 
plane in the (Fe)n/(PbTi03)g/(Pt)8 superlattice. The insets 
illustrate the relative position of the interfacial Ti02 plane 
with respect to the adjacent Fe and Pt atoms. The arrows 
mark the direction of the polarization in PbTiOs. 



DFT-LDA calculations. For two metal electrodes, bee 
Fe and fee Pt, we find lattice parameters of 3.904 and 
3.901 A, respectively, to be compared to the experimen- 
tal values of 4.05 and 3.92 A. We terminate the FEs with 
TiC>2 planes and place the Fe and Pt atoms atop oxygen. 
For each of the superlattices, (Fe)n(PTO)g(Pt)8 and 
(Fe)g(BTO)i3(Pt)io (the subscripts indicate the number 
of atomic layers) , we start with the FE displacements cal- 
culated for the bulk phase and fully relax the structure 
until the maximum force is below 10 meV/A, keeping the 
in-plane lattice constants fixed to that of the correspond- 
ing bulk tetragonal phase of the ferroelectric. In the case 
of BTO, this constraint leads to in-plane tensile strain 
of about 1% and 1.2% in Fe and Pt, respectively, and in 
the case of PTO compressive strain of about —1% and 
—0.9% in Fe and Pt, respectively.— 

We start with discussing the Fe/PbTiOs/Pt superlat- 
tice. The macroscopically averaged electrostatic poten- 
tials of the relaxed supercell for two opposite directions of 
FE polarization are shown in Fig. [51 Within a standard 
FE capacitor model^ the depolarizing field in the ferro- 
electric film is proportional to 8ttP\m /c^fe, where Am 
and dpE are the screening length of the metal electrodes 
and the thickness of the FE film. Within this model, the 
depolarizing field vanishes only if Am — ► 0, in which case 
the surface polarization charges are fully compensated 
by the free charges in the metal. Fig. [2] shows a residual 
depolarizing field of about 0.06 V/A in the PTO layers 
that points in the direction opposite of polarization, in- 
dicating incomplete screening by the metallic electrodes. 
The residual field in the PTO superlattice is about 17% of 
the depolarizing field calculated for an isolated PTO film. 



This ratio is comparable to 2Am jdpTO ~ 13% where the 
FE thickness dpxo is 16 A and screening length of Fe Am 
is about lAJ£ 

Using the Born effective charges— and ionic displace- 
ments relative to the ideal positions, we estimate the 
polarization of the central unit cell to be about 75 and 
70/iC/cm 2 when the polarization points towards Pt and 
Fe, respectively. These values are slightly reduced com- 
pared to the calculated bulk value of 84yuC/cm 2 . (For 
the BaTiOs based structure, we have found 20 and 18 
/xC/cm 2 for the respective polarization, whereas the bulk 
value is 25/LtC/cm 2 .) Fig. [3] shows the rumpling param- 
eter, defined as the relative displacement between the 
cations and oxygen ions for each atomic plane along the 
stacking direction of the PTO superlattice. The ampli- 
tudes of the interfacial rumpling are roughly 0.1 A higher 
when the polarization points towards Pt than when it 
points towards Fe. The higher amplitude in the former 
case is consistent with that observed in Pt/PbTiOs/Pt 21 
and can be most likely attributed to the formation of 
metallic bonds between Pt and Ti at the Pt/PTO inter- 
face and the Fe-0 bonding at the Fe/PTO interface. 

At a FE/metal interface, the ferroelectric polarization 
terminates and induces screening charges in the metal. 
The total screening induced charge density satisfies the 
Poisson equation 

f^ = -(e/e )[6nl(x)+6nt(x)], (1) 

where V c is the Coulomb screening potential in the metal 
and (5n CT is the screening induced charge density of spin 
a =1, J.. If the metal is a FM, then the induced charges 
are spin dependent due to the exchange interactions 
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FIG. 4: The layer by layer magnetic moment changes rela- 
tive to the bulk of two Fe atoms per lateral unit cell of the 
(Fe)ii/(PbTi0 3 )9/(Pt) 8 superlattice. The PTO is in a (a) 
ferroelectric phase where the FE polarization directions are 
marked by the arrows and (b) paraelectric phase. (c)The 
magnetic moment changes calculated within LDA+U (where 
(7 = 8 eV). 

and a local magnetization can be induced at the inter- 
face. The local induced magnetization and the Coulomb 
screening potential at a FM/dielectric- or FM/FE inter- 
face^ can be related through 

5n\x) - 6n±(x) = - ° eV c (x), (2) 

where M and N are the spontaneous magnetization and 
the total density of states, and J is the exchange splitting 
in the FM. Thus the induced magnetization switches sign 
when the orientation of the FE polarization switches. 

To calculate the interfacially induced charge density of 
Fe due to the spin dependent screening, we project the 
density of states (DOS) of the Fe/PbTi0 3 /Pt superlat- 
tice onto each atomic layer of Fe parallel to the interface. 
By integrating the DOS below the Fermi level separately 
for the spin-up and spin-down components, we obtain the 
layer by layer magnetization densities n^{x) — n^(x) in 
the Fe unit cells. In Fig. 0] (a), we plot the change of 
the layer by layer magnetization density with respect to 
the bulk value (~ /lateral unit cell counting two Fe 
atoms) for two opposite polarization directions in PTO. 
The total magnetization of Fe is significantly suppressed 
with respect to the bulk at the Fe/PTO interface and 
slightly enhanced at the Pt/Fe interface. To understand 
this behavior, we plot 6M(x) of Fe when PTO is in a 
paraelectric (PE) phase, i.e., P = 0, in Fig. 4(1)) for com- 
parison. It is clear that similar changes of the magnetic 
moment in Fe at both interfaces exist even in the absence 
of the FE polarization! We start by explaining this effect. 

The reduction of the iron's magnetic moment at the 
Fe/oxide interface is caused by charge transfer from Fe 
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FIG. 5: The projected density of states (PDOS) of interfacial 
Fe atoms in the Pt/PTO/Fe superlattice at the (a) Pt/Fe 
interface, (b) Fe bulk region, and (c) Fe/PTO interface. Panel 
(d) shows the PDOS of the interfacial Ti0 2 layer in PTO. The 
up and down arrows represent the spin majority and minority 
channels, respectively. 

to the oxide caused by several mechanisms. In Fig. [5l 
we show the DOS of Fe at both interfaces and compare 
them to that in the interior bulk region. At the Fe/PTO 
interface as shown in Fig.[5tc), the majority spin state at 
the Fermi level dominates, same as in bulk Fe. Thus los- 
ing charge from Fe to the oxide suppresses the majority 
spin more than the minority and leads to a reduction in 
the magnetic moment of Fe. The outflow of charge from 
Fe into the oxide is due to by two effects. First and the 
larger of the two is the formation of Fe-0 bonds at the 
interface. Electronegative oxygen pulls the charge away 
from Fe. Second effect is filling of the so-called metal 
induced gap states (MIGS) that are clearly seen in Fig. 
EJd). This is universal for any metal/insulator interface 
when the metal Fermi level happens to be in the gap of 
the insulator. At the Fe/Pt interface, on the other hand, 
the density of states of Fe at the Fermi level is reversed 
with respect to the bulk as can be seen in Fig. [5] (a) 
and (b). Thus charge transfer from Fe into Pt due to 
the work function difference (the work function of Fe is 2 
V less than Pt) results in losing more minority spins and 
therefore increases the magnetic moment in Fe. We find a 
similar effect in our Pt/Fe bilayer calculation. The com- 
parison with the P = case corroborates that the large 
moment change relative to the bulk is independent of 
polarization. In the following, we will remove this effect 
from the polarization dependent induced magnetization. 

Before proceeding, we make two relevant comments. 
The charge transfer from Fe to oxide is exacerbated by 
the reduction of the PTO band gap to about 2 eV within 
the LDA, compared to the experimental gap of 3.5 eV. 
By applying a Hubbard U correction to Ti 3d states, the 
band gap of the oxide is opened up22 and the overall 
magnetic moment reduction has been reduced as shown 
in Fig.[3Jc). However, the polarization dependence of the 
induced magnetization which is the main focus of the pa- 
per qualitatively stays the same. Henceforth we only dis- 
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FIG. 6: (a) The magnitude of the polarization induced mag- 
netization for two Fe atoms per lateral unit cell vs. distance 
from the Fe/FE interface for the PTO and BTO based su- 
perlattice. Only half of the Fe cell close to the oxide is plot- 
ted, (b) The electrostatic screening induced magnetization 
SM e i (diamonds) and the magnetization induced by a chemi- 
cal component related to polarization switching SM C (stars). 
The inset shows the induced charge density near the Fe/PTO 
interface. 



cuss the results from the LDA calculations. We further 
note that a recent study^ of a symmetric Fe/BaTiOa/Fe 
junction reported large induced magnetic moments on in- 
terfacial Ti atoms, even higher than the induced moments 
on the interfacial Fe atoms. The authors attributed the 
moment changes on Ti to a significant hybridization be- 
tween the 3d orbitals of Ti and Fe at the interface. In 
contrast to the DOS in Ref.— , our interfacial Ti d or- 
bitals lie mostly above the Fermi level as shown in Fig. 
[5jd). Therefore the moment changes on the interfacial 
Ti atoms are much smaller than those of Fe in our cal- 
culation and will be ignored in our study. 

To focus on the polarization related changes of magne- 
tization we subtract SM(P = 0) from SM(±P) shown in 
Fig. |4j Thus we obtain the polarization induced magne- 
tization of Fe near the PTO interface for two opposite FE 
polarization directions as shown in Fig. [B^a) . We include 
the results from the BTO superlattice as well. In cither 
case, the induced moments decay rapidly as function of 
distance away from the FM/FE interface, and change 
sign as the FE polarization direction switches. This sug- 
gests that if bias is applied to switch the FE polariza- 
tion, the magnetization becomes electrically controllable 
as consistent with our model calculation.— It is interest- 
ing to note that SM, however, does not completely follow 
an exponential function as typically assumed in screen- 
ing models. As the distance from the interface increases, 
the sign of SM oscillates. We attribute this behavior 
to Fridel-like oscillations in the metali^i It is consistent 
with the induced charge density distribution (see the in- 



set) calculated from the electrostatic potential near the 
interface, which also oscillates. 

Unlike the rest of the Fe layers, the induced moment on 
the first interfacial layer of Fe is significantly asymmetric 
with respect to the direction of the polarization. To ex- 
plain the asymmetry, we make a reasonable assumption 
that the total induced magnetization can be decomposed 
into two separate terms, SM e i and SM C . The first term 
is related to electrostatic screening which changes sign 
as the polarization switches and the second term is an 
additional contribution which is related to the chemical 
bonding at the interface, 



6M(+P) = SM el + SM cl 
SM(-P) = -SMei + SM c2 . 



(3) 



Since the opposite polarization directions result in differ- 
ent interfacial atomic configurations (Fe-O-Ti vs. Fe-Ti- 
O as illustrated in Fig.[3]insets), SM C could be different in 
these equations. However, in an attempt to separate the 
effect of the electrostatic screening and chemical bond- 
ing, we assume the difference between SM c i and SM c2 to 
be much less than the screening induced moment change 
and therefore SM C \ ss 5M C 2- We calculate SMp and 5M C 
by subtracting and adding SM(+P) and 8M(—P) and 
plot them in Fig. [BJb) . Note that SM C exists only in the 
first interfacial layer. The magnitudes of both SMp and 
SM C are strongly dependent on the magnitude of the po- 
larization in the FE oxide. In the PTO and BTO super- 
lattice, the maximum magnitude of the screening induced 
moment change SM e i{x) is about 0.4/is and 0.13/is per 
surface unit cell, corresponding to a magnetization den- 
sity of 53 and 14 /iC/cm 2 , respectively. The ratio be- 
tween the two resembles closely the ratio between the 
calculated polarization of PTO and BTO. Thus the elec- 
trostatically induced magnetization is clearly linear in the 
FE polarization, in agreement with Eq. (1) and (2). This 
also suggests that our approximation for 8M C is reason- 
able. On the other hand, the magnitudes of SM C are quite 
different, 0.19 and 0.03/i s , at the Fe/PTO and Fe/BTO 
interfaces, indicating a different interfacial bonding char- 
acter between two different FEs with metals^ 

To estimate the ME coupling, we use the dimension- 
less ratio between the screening induced magnetization 
density and bulk electric polarization of the FE^ i.e., 
T] = SIi p l ■ The maximum values of r\ are estimated to 
be 0.66 and 0.61 at the Fe/PTO and Fe/BTO interface, 
respectively. For comparison, we also quantify the ME 
coupling using the definition that is strickly only suit- 
able for bulk materials HoSM e i = aE Cl where E c is the 
coercive field of the FEs and is the vaccum permibility. 
Using a typical value of E c ~ lOOkV/cm, we estimate the 
maximum interfacial magnetoelectric coefficient a to be 
3 x 10 _10 Gcm 2 /U which is similar in order of magnitude 
as that reported by Duan et al£ However the coupling 
in our study arises solely from electrostatic screening, in 
contrast to the dominant contribution from orbital hy- 
bridization at the interface. 7 
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In summary, we have carried out a first principles 
study of the interfacial magnetoelectric coupling in tri- 
component FM/FE/NM superlattices. Through a de- 
tailed analysis of the magnetic moment changes in Fe we 
have identified the major contributions to the interfacial 
magnetization. First of all, even in the absence of the FE 
polarization, the magnetic moment in Fe is significantly 
reduced by charge transfer from Fe to the oxide owing 
both to the Fe-0 bond formation and metal induced gap 
states in the oxide. The latter effect is somewhat excer- 
bated by the underestimated LDA band gap, that can 
be corrected at e.g., the LDA+U level. Importantly, our 
LDA+U calculations show that once the FE polarization 
is turned on the difference between the magnetic mo- 
ments for two opposite FE polarizations is not sensitive 
to the band gap. The polarization induced change of the 



magnetization at the interface has two separate contri- 
butions. We identify a "chemical" contribution that is 
asymmetric in FE polarization owing to the difference 
in the interfacial atomic geometry, and thus is interface 
specific, and the symmetric contribution arising from the 
spin dependent screening. The last effect is electrically 
controllable, linear in the FE polarization and leads to a 
substantial magnetoelectric coupling. 
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